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INTENTIONALLY  BLANK 


On  the  definition  and  calculation  of  efiective  impedances  and  averaged 
boundary  conditions  with  application  to  the  analysis  of  high  complexity 

frequency  sdective  surfaces 


A  J  Mackay,  J.G. Gallagher 


1.0  Introduction 

Frequency  Selective  Surfaces  (FSS)  have  a  large  number  of  appbcations.  Examples  include 
the  design  and  construction  of  radomes  and  antennas  and  the  construction  of  passive 
microwave  and  millimeter— wave  filters.  With  the  improvements  in  computer  technology, 
numerical  methods  have  become  ever  more  important  in  FSS  design.  However,  it  still 
remains  the  case  that  FSS  with  highly  complex  patterns  are  not  amenable  to  modal 
ar  alysis  methods  in  a  realistic  period  of  time. 

In  the  past,  approximate  methods  have  been  more  frequently  used  than  today,  often 
because  there  was  no  realistic  alternative  available.  Despite,  and  sometimes  because  of, 
modern  computer  abilities,  it  is  well  worth  re-examining  some  of  these  methods  in  the 
light  of  new  requirements.  Indeed,  such  methods  can  often  highlight  aspects  of  the  physics 
that  may  often  be  overlooked. 

In  this  report,  we  focus  our  attention  on  the  use  of  approximate  boundary  conditions 
and  their  use  in  the  analysis  of  complex  structures  within  the  context  of  modal  analysis 
programs  for  FSS  design.  This  may  be  of  particular  importance  in  the  study  of  fractal  and 
fractal-4ike  periodic  patterns  over  restricted  bandwidths. 

Approximate  boundary  conditions  have  been  studied  extensively  by  Kontorovich  et 
al.  (see  [1])  in  the  early  sixties  and  before.  The  motivating  philosophy  assumes  that  it  is 
not  necessary  to  calculate  the  ‘exact’  currents  and  fields  induced  in  the  neighborhood  of  the 
scatterer,  on  a  scale  very  much  less  than  a  wavelength,  in  order  to  determine  far  fields  and 
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transmission  and  reflection  coefGcients  far  &om  the  scatterer.  All  that  is  required  is  to  have 
a  relationship  between  the  averaged  electric  and  magnetic  fields  on  a  scale  possibly  much 
larger  than  the  fine  detail  of  the  scatterer  although  still  small  compared  to  a  wavelength. 
Indeed,  this  is  precisely  what  occurs  with  the  use  of  natural  and  artificial  dielectrics  where 
fields  on  an  atomic  scale  are  not  required  for  macroscopic  field  calculations. 
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1.1  Averaged  Boundary  C!onditions 


Kontorovich  et  al.  [1]  applied  averaged  boundary  conditions  to  the  determinatior  of  the 
reflection  coefficient  of  a  plane  wave  from  a  plane  wire  mesh.  This  was  generalised  further 
in  [2]  and  extended  by  Astrakhan  [3].  In  this  section,  we  define  what  we  mean  by  averaged 
boundary  conditions,  providing  a  generalisation  to  the  definition  given  in  [2],  and  describe 
the  scale  relation  between  averaged  boundary  conditions  and  currents. 

Suppose  we  consider  the  planar  scattering  problem,  appropriate  to  FSS,  in  which  there 
exists  a  complex  and  finely  detailed  structure  in  the  plane  z=0.  The  microstructure  is 
taken  to  be  an  arbitrary  pattern  of  arbitrary  complex  impedance  (generally 
inhomogeneous)  of  zero  thickness  measured  in  ohms  per  square.  Let  us  consider  a  region  of 
surface  small  compared  to  a  wavelength,  but  large  compared  to  the  scale  of  the  detail  of 
the  structure,  such  that  there  is  negligible  phase  variation  across  the  area  of  interest.  This 
can  be  represented  by  an  area  AA,  as  illustrated  below,  with  the  understanding  that  the 
boundaries  are  not  "hard"  but  defined  by  a  weight  function  w^(x)  as  defined  in  the 
following  section. 
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Within  any  part  of  the  structure  outlined  by  A  A,  the  fidds  may  be  determined  in  terms  of 
an  incident  field  by  quasi— static  means.  If  the  structure  is  split  into  distinct  cells,  each  of 
size  and  shape  AA  then  the  i^^  such  cell,  AAi,  is  excited  by  a  field  from  sources 

outside  of  the  cell. 

In  general,  it  is  possible  to  define  the  averaged  fidds  with  respect  to  a  wdghting 
function,  w^(^,  of  the  same  functional  form  throughout  the  z=0  plane.  Ignoring  the  i 
suffix,  we  may  thus  define  the  averaged  dectric  fields  and  currents  over  the  cell  by, 

|‘(X)  =  J^rjwj(x-x')  E'’(x')dx'  (l.la) 

c 

I*(i)  =  i;  (l-lk) 

U 

c  c 

where  and  are  the  electric  fields  and  currents  averaged  over  a  much  finer  scale. 
If  an  exact  point  wise  description  of  the  fields  and  currents  is  defined  and  given  by  E(^ 
and  J(3^,  then  we  may  write, 

U 

=  ^jj w^(x-x')  J(x')dx'  (l.ld) 

U 

where  w^(x)  represents  the  fine  scale  weight  function.  In  general,  however,  E(x)  and  J(x) 
are  not  observables  and  may  not  be  measured  in  any  realisable  experiment. 

In  the  above  definitions  the  domain  of  integration  U  is  over  the  entire  z=0  plane 
and  the  area  functions  A  a  and  Aa  are  defined  by 
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(l.le) 


Aa  =  jj  Wi(x')dx' 

U 

<ya  =  jjw^(x0dx'  (1-lf) 

U 

where  the  weight  functions  w^(x)  and  w ^(x)  are  defined  to  be  local,  real  and  non— negative. 
We  define  such  a  w(x)  according  to  the  rules, 

(1)  w(3^  >  0  Vx€U 

(2)  jj  w(x)dx  exists  and  is  finite 

U 

(3)  w(^  ~  as  |x| 

This  last  condition  indicates  that  w(^  is  chosen  to  be  a  function  that  decays  faster  than 
l/lx|  for  |x|  sufficiently  large.  In  addition,  we  further  restrict  w(3^  such  that 

(4)  w(3^  takes  a  single  maximum  value  when  x=0 

(5)  we  define  the  shape  boundary  as  the  locus  b=b(^),  parameterised  by  6, 
such  that  lw(y/w(0)|=0.5  and  restrict  w(^  so  that  |b|  falls  within  the 
localised  limits  B/2<  |b(^)  |<3B/2  for  some  B>0  V^. 

With  this  last  definition  we  can  specify  that  b^  >  b^,  associated  with  w^(^  and 
respectively,  distinguishing  between  the  two  different  length  scales  involved  in  the 
definitions  of  the  averaged  currents  and  fields.  Furthermore,  we  can  now  associate  Aa  with 
the  area  of  the  region  A  A,  a' d  the  shape  of  A  A  with  the  locus  b(^). 

We  now  seek  to  define  a  general  tensor  impedance  relating  F^(x)  and  at  the 
point  X  by. 


(1.2) 
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with  implied  summation  over  the  index  j.  In  this  planar  problem,  we  only  need  to  consider 
field  and  current  components  within  the  x— y  plane  and  Z=Zij  may  be  represented  as  a  2x2 
matrix.  We  note  that  in  general  Z  will  be  a  function  both  of  the  weight  function  and 
of  the  position  on  the  plane  x. 

Suppose  the  field  incident  on  the  2=0  plane  is  given  by  Eo(^,  then  over  each  AAi  we 
may  assume  Eo(}^  is  approximately  constant.  Let  us  define  the  tangential  component  of 
this  field  over  the  i^i>  cell  as  Ej, .  Then  we  may  write, 

■*  ifndCi)  “  Z(“i(i) .  i)  ihi.)  tor  rreiAi  n  I) 

(1.3) 

where  D  is  that  region  of  the  z=0  plane  where  there  are  non-zero  currents  and  where 
tangential  component  of  the  induced  electric  field  from  currents  everywhere 
in  the  z=0  plane. 


and  G  is  the  familiar  Greens  function  given  by 


e(jr-S')=iS;(k8l-VV) 


^-jkjx-x'l  ■ 

Ij-i'  1 


(1,4) 


(1.5) 


where  ko  is  the  free  space  wave  number,  I  is  the  identity  operator  and  rj  is  the  impedance  of 

free  space.  When  x  6  AAi  but  x  ^  D  (ie  in  the  aperture  regions),  Z(w ^(x) ,  x)  has  infinite 

components  and  J^(x)  has  zero  components  to  match.  If  we  are  considering  the  vector 
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component  of  (1.3)  parallel  to  ,  then  the  right  hand  side  of  (1.3)  is  undefined  in  the 

aperture  region.  However,  since  (1.3)  is  valid  for  arbitrarily  large  but  finite  Z  it  is  useful  to 
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extend  the  region  of  validity  vnth  the  understanding  that  Z(w ^(x) ,  x)  and  (x)  can  not 
be  separated  for  x^D.  Thus  we  now  write 

JJg(x-x')J^(x')  dx'  «  Z(w^(x),^jj(x)  Vx€AAi 
. H  J  tang 


Now,  returning  to  (1.2)  in  combination  with  (1.1a)  and(l.lb),  we  have  the  expression, 


SSi*”®'  +  l;(J “Afi-i')  ifndCiO-ix'  »  J(Vj(x),x).J-jJ Vi(x-x')  I^x')dx' 


Substituting  (1.4)  in  to  (1.7)  gives  us 


E|,fng)  ^  ^JJ w^(x-x')  jj e(x'-x'')J^(x'')  dx''dx' 

■  U  D 

*  Z(w^(x),x).j-Jj w^(x-x')  T^(x')dx'  VxeAAi  (1.8) 


We  now  integrate  (1.6)  over  the  entire  z=0  plane  with  respect  to  the  weight  function 
w^(x-x').  Since  is  assumed  constant  over  AAj,  we  have 


E‘tang)  +  1-JJ w^(x-x'')^ C(x''-x')Z^(x')  dx'dx" 


a  \-jj w^(x-x")  Z(w^(x"),x").J^(x")  dx"  VxeAAi  (1.9) 


We  may  now  compare  (1.9)  and  (1.6)  to  show  the  scale  relation  between  averaged  currents 
and  impedances, 


c  u 

VxeAAi  (1.10) 


with  the  proviso  that  Z(w ^(x) ,  x)  J^(x)  is  defined  everywhere  according  to  the  definition 
in  (1.6). 


1.2  The  determination  of  Z(w^(x),  ^ 

In  order  to  determine  the  averaged  boundary  condition  Z(w^(x) ,  x)  it  is  necessary 
to  know  both  the  local  boundary  condition  Z(w^(x),  x)  and  the  local  current  distribution. 
In  regions  where  there  is  an  aperture  and  zero  current,  the  local  electric  field  distribution 
on  the  aperture  is  required.  These  distributions  can  be  found  by  an  application  of  the 
method  of  moments  using  the  wave  equation  applied  to  Ai  for  each  i  in  isolation. 
Alternatively,  provided  AAj  is  sufficiently  small,  Laplace’s  equation  may  be  used. 
Assuming,  for  the  moment,  that  each  cell  AAi  contains  mj  unknowns  and  there  are  n  such 
cells  (each  of  which  is  different),  the  full  scattering  problem  involving  all  the  cells  can  be 
solved  by  n  matrix  inversions  each  of  size  mi  (l<i<n),  followed  by  a  single  matrix  inversion 
of  size  n.  This  may  be  several  orders  of  magnitude  faster  than  a  full  application  of  the 
method  of  moments  to  the  union  of  all  the  At  resulting  in  one  matrix  inversion  of  size 
E°(mi).  Furthermore,  if  most  cells,  Ai,  are  similar  still  greater  savings  vnll  be  achieved.  We 

i«  i 

should  note,  however,  that  the  correct  determination  of  the  current  and  field  distributions 
over  AAi  requires  a  larger  region  than  AAi  (see  below)  and  thus  the  number  of  unknowns 
will  be  Pioii  for  each  cell  i,  for  some  /?i  possibly  as  large  as  several  hundred  (  ^i=0(o*) 
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where  a  is  defined  below).  Nevertheless,  provided  that  a  solution  of  a  problem  with  S°(mi) 
unknowns  is  more  costly  than  n  solutions  of  problems  with  /Jimi  unknowns,  which  will 
always  be  the  case  for  sufficiently  large  a  problem,  the  method  is  justifiable. 

In  order  to  determine  what  the  local  current  distribution  is  we  assume  that 
can  be  found,  up  to  an  unspecified  constant  multiplier  o,  by  considering  the  solution  to 

e  +  JJ G(x-x')I^'(x')  dx'  =  Z(w^(x),  ^  VxeAA'j 

AA'nD 

(1.11) 

where  VxeAAi, 

and  i=  I 


and  AA<  is  sufficiently  larger  than  and  centred  about  AAj.  This  establishes  a  ‘guard’  region 
where  all  the  fields  in  AA^  induced  by  currents  in  cells  outside  of  AAj  may  be  taken  as 
.ipproximately  uniform  over  AA^.  We  note  that  this  is  valid  only  if  AA<  is  still  significantly 
smaller  than  a  wavelength,  in  which  case  the  Green’s  function  may  be  replaced  by  the 
quasi— static  approximation, 


G(x-x')«5if^VV 


|x-x' 


(1.12) 


If  the  locus  b'(^)  is  the  shape  boundary  of  AA'  and  b(^)  is  the  shape  boundary  of  AA,  the 
radial  distance  between  the  two  boundaries  p(^)=|b'(^)— b(^)|.  If  p  takes  a  minimum 
value  Pain  at  $=0ain  then  we  define  r'=|b'(^«in)|  and  r=|b(^i,in)|  and  the  size  of  the 
guard  region  may  be  determined  by  the  ratio  a=r'/r>l.  This  is  illustrated  in  the  figure 
below. 
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[ 


(>(«)  I* 


A  minimum 


In  theory,  the  minimum  a  should  be  determined  such  that  any  (quasi— static)  dipole 
source  situated  at  the  point  Xq  (with  a  l/|x— Xop  electric  field  dependence)  anywhere  on 
the  boundary  of  A  A'  causes  a  negligible  difference  in  field  strength  over  any  part  of  the 
boundary  of  AA.  Given  this  criterion  and  defining  (  as  the  maximum  permitted  fractional 
difference  in  field  strength  between  the  centre  of  AA  (at  x=0)  and  its  boundary,  then  we 
must  satisfy  the  inequality. 


1  -i  <  i 

(0-1)2  q2  -  o2 


(1.13) 


which,  for  small  gives  rise  to  the  condition  that  a  >  2f(. 

If  we  allow  a  reasonable  choice  of  ^=0.1,  it  is  seen  that  oa!20  which  is  a  figure  that 
seems  rather  high.  In  practice,  it  may  often  be  possible  to  significantly  reduce  a  below  this 
limit.  This  reduction  is  clearly  important  for  any  numerical  solution  of  (1  11)- 
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2.0  Example  of  the  method 


In  order  to  implement  the  method  it  is  first  necessary  to  choose  suitable  weight 
functions.  For  purposes  of  illustration,  we  choose  a  square  pulse  function  of  width  AX  and 
area  AAa=A^, 


Fig  2.1 

We  now  seek  to  apply  the  use  of  approximate  boundary  conditions  to  a  FSS  which  is 
amenable  to  simple  analysis.  In  particular,  we  consider  a  sparse  array  of  electrically  small 
dipoles  on  a  square  lattice.  Part  of  such  an  array  is  illustrated  below. 


Fig  2.2 

In  this  representation  we  choose  a  grid  such  that  each  unit  cell  of  the  FSS  is  divided  into 
nxn  equally  sized  squares,  one  of  which  includes  a  small  dipole  at  its  centre.  The  unit  cell  is 
square  and  of  size  L  such  that  L=nA.  Suppose  that  we  identify  each  such  cell  with  A  A  and 
wish  to  find  Z(w^(^,Xi)  at  each  of  the  i  (l<i<n^)  ceUs,  where  w^(^  is  the  pulse  function 
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associated  with  AA,  defined  by 


n 


x€  AA 
xjf  AA 


(2.1) 


Assuming  that  the  dipoles  are  sparse,  we  may  take  AA'  as  the  unit  cell  or  any  smaller  cell 
since  there  are  no  ciirrent  sources  other  than  at  the  small  dipoles.  In  particular,  we  may 
define  AA'  =  AA  and  for  xcAA. 

For  a  small  thin  dipole  of  length  2h,  radius  a,  situated  at  the  origin  and  aligned  with 
the  y— axis,  the  current  per  unit  length i  [3]  is  given  by. 


where 


and 


ly 


lo 


(coskoy  -  coskph) 
1  -  coskoh 


(2.2a) 


(f  -  foskoh) 
coskoh 


(2.2b) 


for  a<h 


(2.2c) 


and  where  Eo  is  the  component  of  the  incident  electric  field  in  the  y— direction.  In  the 
case  where  kohcl,  this  may  be  further  approximated  to  give  the  small  dipole  result. 


where 


Ir(y)  *  lo 


(2.3) 


If  we  assume  an  electric  field  incident  parallel  to  the  dipole,  we  define  To=I^(0)  and 


t  note  that  this  is  not  the  current  density  per  unit  area. 
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Eo=F  (^-  Using  the  definition  (1.1b)  with  we  can  see  that 


where 


V  x6  AA 
Vx^AA 


h 

lo  =JJi(E)  <i»  =  i/E^  (•>’-?’)  iiy  =  |lol>i  .koh-O 

H  -h 


(2.4) 


for  incidence  perpendicular  to  the  dipole,  this  integral  is  zero.  For  parallel  incidence  the 
integrated  electric  field  over  AA,  <^oes  not  take  an  approximately  constant  value  for 

xeAA.  In  particular,  for  xeAA  when  hcA/2.  Using  the  expression  for  G(x-x') 

given  in  (1.12),  with  the  quasi-static  simplification. 


VV 


3x(y-y')  j  2(y-y')*-x2  * 
[x2*(y-y' )  2]  5^2  [x2+(y-y' )  2j  5^2 


'x'=0 

and  the  definition  (1.1a)  with  1*^(^=E(^  we  can  see  that  if  we  assume  negligable  width  to 
the  dipole, 


(i)i]  (2-5) 

where 


and 


A/2  h  A/2 

=  / dy"  / dy'  (h -y'*)  / dx" 


-A/2  -h 


-A/2 


2(y%y-y')3  -(x"+x)  2 
[(x"+x)2+(y"+y-y')2j  5/2 
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A/2  h  A/2 

ai*(x)  =  fdy"  fdy'(h2-y'2)  (dx" 
-A/2  -h  -A/2 


3(y»'+y_y')(x"+x) 
[(x"+x)»+(y"+y-y')2J 5^2 


Focusing  our  attention  on  the  y  component  of  (2.5),  since  we  do  not  excite  the  dipole  in  the 
X  direction,  we  may  make  use  of  standard  integrals  to  show  that 


:iW)  )  (2-6) 


where 


Ti'"  =  y+A/2 
Ti~  =  y-A/2 
=  x+A/2 
r=2c-A/2 


(2.7a) 

(2.7b) 

(2.7c) 

(2.7d) 


and  where  the  function  32(’7,0  defined  by. 


(2.8) 


( y) 

In  order  to  illustrate  0 1  (^,  in  figures  2.3a,b,c  we  produce  a  contour  plot  of  this  function 

( y) 

for  h=0.1A,  h=0.2A  and  h=0.3A  respectively.  3i  (x,y)  is  plotted  for  -A/2<x<A/2  and 
-A/2<y<A/2  and  A  is  taken  as  unity  (A=l).  It  is  immediately  apparent  that  there  is  a 
significant  variation  in  averaged  field  strength  over  the  area  AA  with  respect  to  the  y 
coordinate.  The  variation  in  the  x  direction  is  much  less  significant.  However,  it  is  clear 
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Figure  2.3c  h=0.3A 
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that  for  hcA,  3i  (^/A2-»o  Vx.  This  is  merely  saying  that  the  averaged  scattered  electric 
field  from  a  small  dipole  becomes  smaller  as  the  dipole  size  is  reduced.  It  is,  for  example, 
relatively  easy  to  show  that, 

^o=i  (j»?o/’rlto)Io  (2h/A)  as  h/A-*0.  (2.9) 

We  may  now  use  (1.10)  to  define  the  terms  of  the  effective  impedance  matrix  Z(w^(^,^ 
for  all  X  within  the  cell  A  A  containing  the  dipole*,  Z  =  z"  ]  ’  ^yy  given  by 

jjEf'da 

^yy(*)  - - z -  Vx€AA 

Javg* 

Assuming  E^(^  is  sufficiently  small,  as  discussed  above  for  small  dipoles,  we  may  use  (2.4) 
and  (2.2)  to  obtain  the  expression, 

Ann  fA'I  2 

Zyy  »  — j  j— jj-  for  h<A/2,  kohcl  (211) 

where  Zyy  is  independent  of  x  for  x€AA.  We  note  once  again  that  this  effective  impedance 
is  approximately  correct  provided  we  define  the  weight  function  w^(^  to  be  the  pulse  basis 
function  given  by  (2.1)  and  thus  take  the  current  distribution  to  be  approximately  uniform. 
Such  an  averaged  current  distribution  coupled  with  the  effective  impedance  must 
automatically  be  consistent  with  Maxwell’s  equations.  This  is  confirmed  in  this  example, 
since  a  small  area  of  material  with  high  self  impedance  (note  that  Zyy  is  large  and 
capacitive  when  hcA  and  kohcl)  is  known  to  exhibit  an  induced  current  distribution  that  is 

*  all  the  terms  of  Z  for  cells  AA  which  do  not  contain  a  dipole  are  infinite. 
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uniform  over  most  of  the  material. 

In  order  to  check  the  validity  of  this  expression,  we  compare  the  far  field  vector 
potential  generated  by  the  induced  current  in  the  dipole  with  the  effective  impedance 
method  applied  to  the  uniform  impedance  square. 

The  scattered  far  field  generated  by  an  electrically  small  dipole  with  negligible 
width  is  defined  by  a  vector  potential  Adip(x).  Using  the  expression  in  Jackson  [4],  the  far 
field  vector  potential  is  given  by 


Lim  A  W  = 

r^m  "dip  ^ 


-h 


where  Jdip(20  is  the  current  distribution  on  the  dipole.  This  integral  has  already  been 
evaluated  in  (2.4),  thus 


Lim  A  (x)  = 
-dV-" 


e-jkr 

47r 


4 

3 


Ioh£  = 


e-j'^^  1  jEg 
ITT  "15  F 


(2.12) 


Similarly,  the  scattered  far  field  from  the  small  impedance  patch  with  constant  surface 


impedance  Zyy  and  pulse  current  distribution,  Jeff=Jo£  ,  gives  rise  to  a  far-field  vector 


potential. 


i.e. 


A  A 
5  5 


Lim  A  (jO  = 


JJdxdy 


-A  -A 

7  7 


Lim  A  (^  =  Jo£ 

r-tQD  ?TT 


(2.13) 
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but  the  impedance  boundary  condition  states  that 


Jo  8  Eg^VZyy 

assuming  that  the  scattered  electric  field  over  the  small  impedance  square  is  small 
compared  with  the  incident  field  strength  Substituting  the  expression  (2.11)  into 
(2.13)  shows  that 

LimA  =  koh=»  (2.14) 

r-4CB 

which  is  the  same  as  (2.12),  confirming  the  comparison.  We  would  also  note  that  in  this 
example  we  assume  that  the  scattered  field  is  much  smaller  than  the  incident  field  and  thus 
an  FSS  composed  of  such  dipoles  or  impedance  squares  suffers  negligible  mutual 
interactions.  Therefore,  we  conclude  that  the  properties  of  the  FSS  arrays  of  dipoles  and 
the  effective  impedance  squares  must  also  be  identical. 


17 


3.0  Conclusion 


Although  we  have  not  presented  an  example  showing  a  use  of  the  theory  of  averaged 
boundary  conditions  for  a  useful  FSS,  we  have  illustrated  that  the  method  is  valid  and  have 
developed  a  generalisation  of  the  theory  over  existing  work  [1—3].  It  is  hoped  that  the  work 
presented  here  provides  a  basis  for  the  development  of  predictive  software  tools  to  allow 
the  analysis  of  structures  which  are  currently  too  complex  to  be  examined  using 
conventional  modal  analysis  methods. 
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